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We study the competition between phase definition and 
quantum phase fluctuations in interference experiments be- 
tween independently formed Bose condensates. While phase- 
sensitive detection of atoms makes the phase progressively 
better defined, interactions tend to randomize it faster as the 
uncertainty in the relative particle number grows. A steady 
state is reached when the two effects cancel each other. Then 
the phase resolution saturates to a value that grows with the 
ratio between the interaction strength and the atom detection 
rate, and the average phase and number begin to fluctuate 
classically. We discuss how our study applies to both recently 
performed and possible future experiments. 
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Recent macroscopic interference experiments |lj|2|] have 
revived interest on the question of whether two inde- 
pendent Bose condensates can have a well-defined phase 
HQ. Guided by the experimental result that a precise 
phase was indeed observed, the ensuing theoretical work 
has noted that, even if the initial relative phase is 
random, as happens for pairs of condensates prepared 
in a Fock state, the phase becomes progressively well- 
defined as atoms are emitted from the compound system 
and recorded in a phase-sensitive detector. The resulting 
view accords well with the standard quantum measure- 
ment picture according to which the phase is created by 
the very act of measurement; if one tries to measure the 
phase one indeed observes a definite phase. However, 
these works have not discussed how the phase definition 
process is affected by the mean-field atomic interactions, 
which play a competing role by causing quantum spread- 
ing of the phase. In this article we study the role of in- 
teractions, and show that they modify the above picture 
considerably. While the final conclusion remains true 
that a well-defined phase is expected in a wide range of 
experiments, we argue that this is possible only because 
the interactions are sufficiently weak, and calculate ex- 
plicitly how weak they must be in a specific experimental 
setup. The central idea is that, due to interactions, phase 
is not a constant of motion. An experiment will be ef- 
fective in measuring the phase only if the measurement 
is sufficiently intense (in a sense that will be made pre- 
cise) to overcome the phase internal dynamics opposing 
definition. 

We study the effect of interactions within the frame- 
work of the interference experiment considered by Castin 
and Dalibard [Hi. Two atom beams emitted from inde- 



pendent atom lasers || are made to converge on an atom 
beam splitter |lC| and the emerging beams are directed to 
two detectors. This experiment has the attractive feature 
that, without being unrealistic, it lends itself to a funda- 
mental theoretical analysis, since it captures the essence 
of phase measurement in its most simple form. With a 
good understanding of the central physics, we will see 
that it is possible to draw some semiquantitative conclu- 
sions on interference experiments, such as that already 
performed at the MIT §. 

We consider two independent Bose condensates, each 
confined in a harmonic trap. A Hamiltonian for the mean 
field description of the interaction is easily derived from 
a Taylor expansion of the energy of each condensate [[ll] . 
We have N a (Nb) atoms in condensate A (B), N — N a + 
N b and n = N a - N b . For \n\ < N, 



H = 2E(N/2) + E c n 2 /2 



(1) 



describes the coherent dynamics of the total system. 
E c = E"(N/2) determines to lowest order in n the self- 
interaction within each condensate, and can be derived 
from the ground state properties of an interacting Bose 
gas [ ft2"|| . In a quantized version, the atom numbers N a 
and Nb are given by the operators a' a and b^b. 

Following Ref. §, we consider an interference exper- 
iment in which the condensate atoms are emitted from 
the traps and guided to a 50-50 beam splitter. The ac- 
tion of the two detectors D± on the system is described 
by the operators 



(2) 



where 7 is the outcoupling rate per atom, and the expec- 
tation values I±(t) = (ip(ti\Cl.C±\Tp(t)) give the corre- 
sponding detection rates [^3| . 

A convenient representation is provided by the phase 
states 
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If the wave function of a state \ip) in the phase representa- 
tion ip{4i) = jsi{4>\ip) exhibits a sharp peak at an average 
value 4>, it is easy to see that the detection rates read 
I± = lN[\ ± cos(2</>)]/2. This property of full fringe visi- 
bility makes the phase states (p|) ideal for interferometric 



purposes. The number difference operator n 
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acts on the phase states as n|0)jv = — i(d / d<f>)\(j>) ^ . 
Thus n and (f> are conjugate quantities, and, for a sharp 
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phase distribution, the inequality AnAtp > 1/2 applies. 
The Heisenberg equations of motion d<p/dt = —E c n/Ti, 
dn/dt — 0, lead to a time-dependent interference pattern 
I± (t) oc 1 ± cos(2E c ht/h). 

A necessary feature of phase-sensitive detection is that 
it does not give information on which condensate the 
recorded atom comes from. Each detected atom may 
come with equal probability from either condensate A or 
B. As a result, the uncertainty in the relative atom num- 
ber after k detections grows like An = Vk. In the ab- 
sence of interactions Q , a minimal wave packet is formed, 
so that A(f> = \j2\fk. Let us assume that the same result 
holds true in the presence of interactions. Then we can 



conclude that the increase An 2 — > Ar 



1 in the num- 



ber uncertainty occurring in each additional detection, is 
accompanied by a corresponding decrease in the phase 
uncertainty, 
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On the other hand, interactions generate phase dynam- 
ics. The mean value of the phase shows a drift, cf>(t) = 
<f>(0) — E c nt/h, while its variance undergoes ballistic 
spreading, A4>(t) = A0(O) + E c Ant/h. Thus between 
detections the phase uncertainty grows at a rate propor- 
tional to An, which in turn grows with the number of de- 
tections. Since is the total detection rate, the mean 
increase of phase variance between detections is 



Acf> — * A4> ■ 
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(5) 



In the early stages of the measurement process An is 
small and the effect of interactions is negligible; the vari- 
ation in ([5]) is much smaller than that in (^). However, 
as k grows, the two effects become comparable and can- 
cel each other. Then a steady state is formed in which 
the phase becomes defined by the detections at the same 
rate at which it is randomized by the interactions. The 
stationary values of the variances are 



^ = (l) V4 > A 
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Saturation occurs after k* — An 2 atoms have been de- 
tected. In the absence of interactions, Atfi s — 0, and the 
phase can become arbitrarily well-defined. In contrast 
to that, Eq. (g) indicates that interactions cause an in- 
trinsic limitation in the achievable phase resolution. A 
steady state with a well-defined phase A<f> <C 2tt can 
be established only if interactions are sufficiently weak 
(k <C 1). For a 3d harmonic trap in the Thomas- Fermi 
limit, E c ~ 7V-3/5 [JyjrflUlJ] an d thus the effective inter- 
action parameter, n ~ iV -8 / 5 , increases during the de- 
tection process. However, we shall focus on time scales 
where N is practically constant (k <C N). 



For a more quantitative analysis we have performed 
a numerical simulation of the time evolution in a sin- 
gle experiment based on the quantum jump method p3[ . 
Starting from an inital Fock state |^(0)) = \N a , Nt>), the 
time evolution over a sufficiently small time interval At 
can be simulated as follows. If an atom is detected at 
D±, the wave function changes as 



W*)>-c±W*)>- 



(7) 



These detections occur with probabilities 

P±(t,At) = {rp(t)\clC ± \rp(t))^ « 1- (8) 

If no detection takes place, the system propagates from 
t to t + At with the effective Hamiltonian 



H - z(C|C + + ClC_)/2 = H- i-y((Ja + b%)/2. 



Since the propagation with H e g as well as the detection 
(0) is non-unitary, the wave function has to be normal- 
ized after each time step. In the present case, an impor- 
tant simplification is brought about by the fact that the 
time evolution with H e g leaves the detection probabili- 
ties unchanged. This means that the detection times are 
uncorrelated, and the time t between consecutive detec- 
tions is exponentially distributed, w(t) = jN exp(— jNt) 
being its probability distribution. Therefore the stochas- 
tic time evolution described above is equivalent to one 
which can be implemented more efficiently, namely that 
the system propagates with the Hamiltonian (|l|) over an 
exponentially distributed random time t and after each 
propagation a detection (^) is performed at either D + or 
L>_ with probability P ± (t) = I±/(I+ + 2_). 

Figure |l] shows how the number variance evolves 
in a typical experimental run for different interaction 
strengths. Starting the stochastic time evolution with 
the initial state \N a ,Nb), the variance An grows with 
the square root of the number of detected atoms as in 
the interaction free case ||. After k* atoms have been 
detected, i.e. after a time t* = (ANjEc/h)^ 1 / 2 , the in- 
teraction becomes effective, and the variances saturate 
at the values (|^). Figure || shows the obtained An s and 
A(j) s as a function of the interaction strength. For weak 
interactions, our numerical studies confirm to a good ap- 
proximation the estimates @. We also find good agree- 
ment with our asumption of a gaussian shape for the 
wavefunction at large k (see Fig. ||). 

Figure |^ also reveals that for strong interactions (k ^> 
1) An s saturates as a function of k. In this case, we 
expect the phase to become completely random. This 
corresponds to the limit in which interactions dominate 
over phase-sensitive detection. Figure ||a shows, indeed, 
that in this case, instead of a gaussian phase distri- 
bution we find asymptotically a completely undefined 
phase. On the other hand, the random variation of n 
by ±1 makes each detection act like a random kick, in 
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analogy to the kicked rotor problem pl| with n play- 
ing the role of the angular momentum. This dynamical 
system is known to have a regime where the angle dis- 
tribution spreads over the circle at long times, while the 
angular momentum distribution becomes exponentially 
localized within a length that depends only on the kick 
strength, and not on the parameters of the free propaga- 
tion Jl5[ . Thus we expect to find wave functions of the 
form tp n ~ exp(— \n — n|/4An s ), where An s is indepen- 
dent of k. This is confirmed in Figs. |^ and ||b. 

In the non-interacting case |5j each detection amounts 
to a partial measurement of the phase, the measure- 
ment becoming more precise as more atoms are recorded. 
However, what is actually measured is cos(20) through 
the observed intensities I±. Thus the whole measure- 
ment process determines the phase only up to a sign 
(— 7r/2 < <fi < 7r/2), and the condensate is in a super- 
position of two minimal wave packets with peaks at ±</>. 
Within each peak An = \fk and Atfi = (2An) _1 . 

To study the effect of interactions on the phase defini- 
tion process, let us assume first that the initial number 
difference is no ~ VN 3> An s . Figure ||a shows that, in 
contrast to the interaction free case, after an initial tran- 
sient the system evolves into a single wave packet with 
an average phase <f> = E c nt/h. The reason is the follow- 
ing. Since n ~ no, the sign of n is known. On the other 
hand, the detection record yields both, I±(t) ~ cos(2</>) 
and I±(t) ~ nsin(2(/>). This makes the measurement of 
4> unambigous p^ , p7[ |. 

When n < An s , the sign of n cannot be measured 
with enough precision and, as in the non-interacting case, 
4> is determined from the measurement only up to a sign. 
Accordingly, the phase distribution evolves into two wave 
packets. The sign of n is correlated with <j> in such a way 
that the sign of nsin(20) is the same for both compo- 
nents. This can be seen in Fig. ^b, where the two wave 
packets, each with an opposite mean phase value, drift 
with opposite velocity. 

Inspection of Fig. ||a seems to suggest that n stays 
steady once A<p and An reach saturation. However, a 
closer look reveals that this is not the case. We find 
that n remains almost constant during the initial interval 
k < k* (not shown). Once An saturates to An s for k > 
k* , n begins to fluctuate classically with the spread over 
many runs growing like \fk. This implies that the total 
number uncertainty obtained by averaging the quantum 
averages n over many runs grows steadily like yk. This 
should be expected, since the result of such a combined 
average must be identical to the full quantum average of 
two condensates emitting at random. 

To understand the MIT experiment [Q, we may iden- 
tify the emitted atoms with those in the region where 
the two wave packets overlap. Then we may replace the 
outcoupling rate per atom 7 by the rate of overlap in- 
crease, which is of order 10s _1 . Since E c /h w 0.01s -1 
for N w 5 x 10 6 sodium atoms, this results in a phase 



resolution A(f> s «3x 10~ 3 , which is consistent with the 
observation of a clear interference pattern Qj. 

In summary, we have studied the influence of self- 
interaction on the phase definition in a macroscopic inter- 
ference experiment between separate Bose condensates. 
Our main finding is that the ratio between interaction 
strength and detection rate limits the phase resolution. 
This arises when the phase randomizes due to interac- 
tions at the same rate at which it becomes defined by 
the successive partial measurements performed at each 
atom detection. We have derived an analytical expres- 
sion for the phase resolution which has been confirmed by 
a Monte Carlo simulation. Once the interference pattern 
reaches its intrinsic resolution limit, the average phase 
and number fluctuate classically. The studied model al- 
lows us to understand in simple terms why, despite inter- 
actions, good phase resolution can be obtained in macro- 
scopic interference between Bose condensates. It will be 
most desirable to perform experiments that explore the 
fading of the interference pattern with increasing inter- 
action strength. 

We wish to thank C. Bruder and W. Ketterle for 
valuable discussions. This work has been supported by 
the EU TMR Programme under Contract No. FMRX- 
CT96-0042 and by the Direction General de Investigation 
Cienti'fica y Tecnica under Grant No. PB96-0080-C02. 



[1] M. R. Andrews et ai, Science 275, 637 (1997). 
[2] B. P. Anderson and M. A. Kasevich, Science 282, 1686 
(1998). 

[3] P. W. Anderson, Basic Notions of Condensed Matter 

Physics (Benjamin, Menlo Park, 1984). 
[4] A. J. Leggett and F. Sols, Found. Phys. 21, 353 (1991). 
[5] Y. Castin and J. Dalibard, Phys. Rev. A 55, 4330 (1997). 
[6] J. Javanainen and S. M. Yoo, Phys. Rev. Lett. 76, 161 

(1996). 

[7] J. I. Cirac, C. W. Gardiner, M. Naraschewski, and P. 

Zoller, Phys. Rev. A 54, R3714 (1996). 
[8] M. Naraschewski, H. Wallis, A. Schenzle, J. I. Cirac, and 

P. Zoller, Phys. Rev. A 54, 2185 (1996). 
[9] I. Bloch, T. W. Hansen, and T. Esslinger, Phys. Rev. 

Lett. 82, 3008 (1999). 
[10] C. S. Adams, M. Sigel, and J. Mlynek, Phys. Rep. 240, 

143 (1994). 

[11] I. Zapata, F. Sols, and A. J. Leggett, Phys. Rev. A 57, 
R28 (1998). 

[12] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[13] K. M0lmer and Y. Castin, Quantum Semiclass. Opt. 8, 

49 (1996), and references therein. 
[14] G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6 (1996). 
[15] G. Casati, B. V. Chirikov, F. M. Izrailev, and J. Ford, in 

Stochastic Behavior in Classical an Quantum Hamilto- 



3 



[16] 



[17] 



nian Systems, Vol. 93 of Lecture Notes in Physics, edited 
by G. Casati and J. Ford (Springer, Berlin, 1979), p. 334. 
An interesting corollary is that a set of sequential partial 
measurements of the phase <j> yields a partial measure- 
ment of its time derivative <f>, and hence of its canonical 
conjugate n. The compromise is the formation of a gener- 
alized coherent state in phase and number. This conclu- 
sion is consistent with the ubiquity in nature of coherent 
states behaving like classical entities. We note that these 
coherent states do not require the presence of a harmonic 
confining potential as assumed in Ref. |l7[ . 
W. H. Zurek, S. Habib, and J. P. Paz, Phys. Rev. Lett. 
70, 1187 (1993). 



25 



S 
< 



10 

5 




1 


1 1 1 




K = io- 5 ~_ 




K " 10 ' 


I 


Li 

" K = 10" 3 
1 1 1 



100 



200 300 

k 



400 500 



FIG. 1. Time evolution of the number variance An in a 
typical quantum trajectory for different interaction strengths 
k = Ec/h-yN. The full line depicts the behavior in absence of 
interaction. Initially, the total atom number is N = 10 6 and 
the number difference is no = 200. 
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FIG. 2. Stationary variances An s and A(j> s for different 
interaction strength. The average values in the steady state 
(dots) have been calculated from 20 simulations with initially 
N = 10 6 atoms. They are in good agreement with the analyt- 
ical estimates in Eq. ^ (full lines). We do not plot A(f> s for 
large values of k because it is not well defined once it becomes 
of order unity. 
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FIG. 3. Phase (a) and number (b) distribution after 500 
detections for weak (full lines, no = 200) and strong (dashed, 
no = 0) interaction for iV = 10 6 . 




FIG. 4. Time evolution of the phase distribution |?/>(<^))| 2 
for k — 2x 10 -4 and JV = 10 5 . The initial number differences 
are no = 200 (a) and no = 1 (b). For graphical reasons, we 
plot in panel (b) only after every tenth detection. 
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